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SUMMARY 


A method is proposed for the electronic recording of pulsed- laser holograms. By 
controlling the polarization sensitivity of each resolution element of the detector, an 
arbitrary phase can be added to the phase of the image wave . This feat is accom- 
plished by selecting a polarization sensitivity sin [a(x,y)/2] in the x-direction and a 

O 

polarization sensitivity cos [a(x,y)/2] in the y-direction, where a(x,y) is the phase. 
The reference wave is circularly polarized, and the object wave is linearly polarized 
at equal angles relative to the x- and y-directions . The detector* s ability to resolve 
the object wave is not affected by adding a phase in the manner specified. 

If the hologram is scanned with an x-directed scan velocity V and the phase is 
chosen to be given by 


a(x,y) = <5 sin 27rfx 


then a phase 


a(t) = 6 sin 27rfV x t 

appears in the signal generated from the electronically recorded hologram. If the fre- 
quency f is also the maximum spatial frequency in the object wave, the various terms 
of the hologram yield signals in nonoverlapping frequency bands . The signal corre- 
sponding to the self-interference terms of the hologram can be eliminated by bandpass 
filtering. The hologram constructed from the filtered signal yields nonoverlapping 
image waves even though the original hologram was constructed with coaxial reference 
and object waves. 

The pulsed-laser system is based on a similar system for electronic recording 
and processing of continuous-wave holograms. The principle of the continuous-wave 
holography system for electronic recording of holograms and image processing is re- 
viewed. 


INTRODUCTION 

A method has been proposed for the electronic recording and remote transmission 
of holograms (refs. 1 and 2). This method has the advantage that some common pro- 



cedures of coherent optical processing such as cross correlation and spatial filtering 
are readily incorporated in the recording process . 

This paper defines a system that could be used to duplicate the properties of the 
continuous-wave-holography system proposed in references 1 and 2 when the illumina- 
tion is provided by a pulsed laser. Since the heterodyne recording essential to the 
concept is introduced in the spatial domain rather than in the time domain, the method 
is called simulated heterodyne recording. The concept is tested mathematically, and 
the required hologram detector (television camera tube) is defined. 

The systems discussed in references 1 and 2 and in this paper can be used wher- 
ever three-dimensional image processing is required. However, the development of 
a pulsed- laser system is motivated by demands for optical instrumentation for use in 
turbomachinery research and development. The pulsed- laser technique is particularly 
applicable to the holographic flow-visualization methods of interest at the Lewis Re- 
search Center (refs. 3 and 4). For this reason, examples of importance in flow visu- 
alization are discussed. 

This paper is both a review and extension of the work reported in an NASA- 
sponsored thesis (ref. 1). Before the pulsed- laser theory is developed, the material 
in reference 1 is reviewed by discussing the limits and restrictions on electronic re- 
cording of holograms. The analysis of an object wave for its essential information 
content and a method for electronic recording of holograms and image processing are 
presented. The paper begins with a brief section giving the reasons for recording a 
hologram electronically. 


REASONS FOR RECORDING HOLOGRAMS ELECTRONICALLY 

A primary reason for recording a hologram electronically is to permit on-line 
processing and analysis of an optical image. Although holography provides a "picture 
worth a thousand words, " definitive, quantitative results in flow-field studies, for ex- 
ample, usually demand other coherent optical techniques such as laser Doppler veloc- 
imetry. However, with image processing and analysis, specific attributes of an image 
can be evaluated. For example, in principle, the electronic recording method de- 
scribed in this paper could be combined with a flow- visualization system to measure 
shock-wave orientation . 

Another reason for recording a hologram electronically is to display a changing 
three-dimensional scene in quasi- real time . For example, three-dimensional strobo- 
scopic imaging of an intrablade flow field during flutter in a compressor is conceivable. 

Finally, of less importance in the present context is the possibility of combining 
electronic holography and stereoscopy for a "three-dimensional" television display. 
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The multiple-view holography system described in reference 2 would serve as a com- 
ponent of such a display . 

Regardless of the reasons and methods for recording a hologram electronically, 
the low spatial resolution of electronic detectors (television camera tubes) is the pri- 
mary constraint. That constraint is discussed in the next section. 


LIMITS AND RESTRICTIONS ON ELECTRONIC RECORDING OF HOLOGRAMS 

At best, electronic recording of holograms does not have sufficient resolution to 
allow the methods that have made dual-beam holography convenient and routine to be 
used. These methods have been described extensively in the literature of holography 
(e.g., refs. 5 and 6). Because of the completeness of the literature, only a brief dis- 
cussion of dual-beam holography is presented here to show where the resolution re- 
striction asserts itself. 

The dual-beam hologram was proposed by Gabor in 1948 (ref. 7). In recording 
the hologram, two monochromatic beams of light are incident on a plane called the re- 
cording plane, as shown in figure 1(a). The first beam, called the object wave, is 
represented by the complex phasor U(x, y, z). The second beam, called the reference 
wave U R (x, y, z), has the same polarization as the object wave and is coherent with it. 
The two beams interfere at the hologram-recording plane to yield an intensity that is 
proportional to 

I(x, y, z) = | U R (x, y, z) | 2 + | U(x, y, z) | 2 + U R (x, y, z)U(x, y, z) + U R (x, y, z)U*(x, y, z) (1) 

where the asterisk denotes complex conjugation. A so-called linear recording is made 
such that the amplitude transmittance T(x,y, z) of the recording is directly propor- 
tional to the intensity 


T(x,y,z) = yl(x,y,z) (2) 

The reconstruction procedure shown in figure 1(b) consists of illuminating the 
recording with a replica of the reference wave U R (x, y, z). Four waves are recon- 
structed and constitute the four diffraction orders of a linearly recorded dual-beam 
hologram. At the hologram-recording plane, these waves, in order of the terms in 
equation (1), are 


y|U R (x,y, z)| 2 U R (x,y, z) 
T I U (x, y , z )| 2 U R (x,y, z) 
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r|U R (x,y, z)| 2 U(x,y, z) 
yU R (x,y, z )U R (x,y, z)U*(x,y, z) 


Only the third wave is proportional to the original object wave. The first wave is pro- 
portional to the original reference wave . The second wave arises from a term where 
the object wave has interfered with itself. The fourth wave, if the object wave origi- 
nated from a three-dimensional object, forms a real image. 

The resolution requirement of the recording process is determined from the 
Fourier transformation of the third term in equation (1). At the hologram- recording 
plane, where z is constant, the spatial transforms of the reference and object waves 
are defined by 


f f °° — j27rf x-j2irfv 

0(f v ,f v ) = I U(x,y)e r dxdy 




x’ y' 




x’-y' 


( 3 ) 


Both U and U R are assumed to be essentially nonzero over the hologram segment 
only. All interpretations are simple if U and U R are imagined to be composed of 
plane waves 


j27Tf X+j27Tf v 

°( f x’ f y> e ^ 


where 


R(f x »f y )e 


j27rf x x+j27rf y y 


2 7rf = k cos cp 

A A 

= k cos cp 

J 

where k = 2n/X is the wave number and cp and cp are the angles between the 

x y 

propagation vectors and the x- and y-axes, respectively. 
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The transform ^{Ig} of the third term 1^ in equation (1) is given by 


^{Ig} =°(f x ,f y ) ® R*(f x ,f y ) = 


IL 


°< £ X - £ x. K - £ v> B * 


w 


where the symbol (x) denotes the convolution operation defined in equation (4). If O 
and R are approximately band limited, the following conclusions are easily drawn by 
considering the limiting plane waves. Only the f variation is considered since the 
same conclusions hold for the f y variation . The spectra O and R are considered 
to be symmetrically distributed about f QX and f rx and to have half -widths Af ox and 
£f rx , respectively. Then is distributed in the interval 

Tf - f - (Af + Af ), f - f + (Af + Af ft 
Lox rx 1 ox rx ;> ox rx ’■ ox rx'J 

Clearly, the maximum frequency that must be recorded is minimized at Af + 

OX 

Af_„ if f_ v = f_„. Then, the reference and object beams are coaxial. If the beams 

X A UA IA 

are not coaxial, the maximum frequency that must be recorded is increased by the dif- 
ference 


I f ox - f rxl = “ 1 cos ^ox - cos ^rxl ( 5 > 

The virtual-image axis is determined by 

27rf ox = k cos ^ox ( 6 > 

and the real-image axis is determined by 

2,r < 2f rx - W = k cos '"rx < 7 > 

where is the angle of the real- image axis. The axes of the real- image and 

virtual-image waves are coaxial if f „ = f . 

0 rx ox 

Finally, the first two terms of equation (1) can be transformed. The transforma- 
tion of the first term is distributed in the interval (-2 Af rx , 2 Af rx ) and the transfor- 
mation of the second term is distributed in the interval (-2 Af^, 2 Af ). Generally, 
Af-... > Af , so that upon reconstruction it can be shown that self-interference terms 

UA XA 

produce radiation in the spatial-frequency interval 

If - 2 Af - Af , f + 2 Af + Af "T 
|_ rx ox rx’ rx ox rxj 
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But the term yields radiation in the interval 

ff - Af - 2 Af , f + AE +2Af 1 
Lox ox rx* ox ox rxj 

If f = f , it is clear that the self-interference terms yield a larger beam divergence 
rx ox 

than do the image-wave terms . If f = f and Af is small, the maximum fre- 

ox rx rx 

quency that must be resolved is determined by the half -bandwidth of the object-wave 
transform. 

We now have the information necessary to define the resolution required by any 
hologram-recording process. The use of coaxial object and reference waves minimizes 
the maximum frequency that must be recorded and, thereby, minimizes the resolution 
requirement of the recording medium. Unfortunately, then the diffraction orders over- 
lap. The signal-to -noise ratio of the reconstruction is reduced enormously: Only one 
order contains the original object wave, and three of the orders are simply noise. 

At the expense of a greatly increased resolution requirement, Leith and Upatnieks 
proposed their famous solution of this overlap problem (ref. 8). If the object wave is 
distributed about f = 0 (beam axis perpendicular to recording plane), choosing the 

reference -beam angle q> such that 

rx 


f — 3 Af 
rx ox 


( 8 ) 


assures separation of the orders. The Fourier transform of equation (1) is shown 
schematically in figure 2(a) for such a choice. Clearly, the object -wave term is in an 
entirely different spatial band than are the other terms . Most practical applications of 
holography employ the off-axis reference-beam method of Leith and Upatnieks . This 
method makes holography convenient and routine . 

The extreme difficulty of using the Leith- Upatnieks method with electronic record- 
ing is readily shown by an example. Consider a 525-line television camera tube with 
an 18-millimeter diagonal. The tube conveniently has the same number of resolution 
elements in the horizontal and vertical directions. For a 3:4 aspect ratio, the horizon- 
tal resolution f in sinusoidal cycles per millimeter falls in the ordinary photographic 
resolution range, about 18 per millimeter. The reference frequency must be centered 
at 3 Af (fig- 2(a)) and the camera tube must resolve 

OX 

4 Af ox = £ (9) 


Hence, the maximum recordable object -wave frequency Af QX is about 4.5 per milli- 
meter, a value probably not worth the effort. 
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If the hologram is recorded electronically, there is approximately a one-to-one 
correspondence between horizontal spatial frequency and temporal frequency. With 
filtering, it is possible to use a single sideband of an image term to construct an ade- 
quate hologram for a symmetrical object wave. As shown in figure 2(b), the reference 
wave must be centered at 2 Af and the hologram need resolve only 3 Af . The 
maximum detail is now 33 percent of the available resolution, or 6 per millimeter, a 
value still not worth the effort. Clearly, another solution of the overlap-of -orders 
problem must be found if electronic recording is to be used. 

In addition to not having enough resolution for the off-axis reference-beam method, 
an electronic detector may not have enough resolution for all object-wave half-widths 
AL A multiple- view holography method such as the method of reference 2 is a solu- 
tion of this problem. However, a large object-wave bandwidth is usually a matter of 
convenience rather than necessity. In fact, in flow- visualization applications, the 
bandwidth can be set fairly arbitrarily. The information content and bandwidth of an 
object wave are discussed in the next section. 


INFORMATION CONTENT AND BANDWIDTH OF AN OBJECT WAVE 

The principles that govern how well a holographic system records and replicates 
an object wave are the principles of optical imaging discussed extensively by Born and 
Wolf (ref. 9) . The modeling of an object wave such that the essential information can 
be isolated is discussed in reference 1 . Hence, only the results are stated here and 
supported by some examples. The object-wave-processing ability of a holographic 
system can be summarized in a few paragraphs. 

Electronic holography is quite capable of recording the macroscopic patterns of 
surface and volume that identify an object. These patterns have transforms within the 
resolution range of the detector. Electronic holography is also capable of defining an 
object wave in three dimensions. This property is determined by the depth of focus in 
the image space and the depth of field in the object space of an optical system. These 
depths are determined by the numerical aperture of the recording method. In elec- 
tronic holography, the numerical aperture is the sine of the half-angle of the largest 
resolvable pencil of rays from a point source . 

The adequacy of an electronic holography system is really challenged by physio- 
logical rather than physical requirements. Wide-angle viewing of an extended object 
requires a large object-wave bandwidth. Stereoscopic viewing requires that each point 
on an object direct light to both eyes. These requirements are met by a diffuser or a 
diffusing surface whose micro structure increases the object-wave bandwidth. These 
statements are amplified in the following paragraphs . 
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As shown in reference 1, the transform of the object wave from a surface can be 
expressed by the equation 


_ -j7rx(z-d)(r+r ) 

°( f x , fy) = ,y',z')} <g> :f{T(x*,y',z')}Je v x y/ (10) 

where d is the distance from the recording plane to a local origin of the surface, 
I(x',y', z') is an equivalent illumination function, and T(x',y' , z') is an equivalent 
transmittance or reflectance of the surface. The function T (x' , y ' , z’) contains the 
macroscopic variations in phase and magnitude that define the surface pattern. The 
maximum spatial frequency of the object wave is obtained by adding the maximum fre- 
quency contained in the illumination I(x' , y' , z' ) to the maximum frequency contained in 
the equivalent transmittance T (x' , y ' , z' ) . 

In most cases, T(x',y', z') will not have enough bandwidth to use the full numerical 
aperture of an electronic detector. Then a diffuser must be employed to increase the 
bandwidth of I(x' ,y', z') for minimum depth of field. In fact, a diffuser should be used 
with a bandwidth determined by the full numerical aperture of the detector. Because 
sum and difference frequencies are involved, the information in T(x',y',z') is still re- 
tained. As explained by Leith and Upatnieks (ref. 10), the required dynamic range of 
the recording process is reduced by using a diffuser. Liu and Gallagher (ref. 11) ex- 
plain how to design an optimum diffuser for an object. 

These comments are reinforced by a simple, yet relevant, example diagrammed 

in figure 3. A discontinuity in refractive index (a shock wave) is illuminated with a 

plane wave. A hologram is recorded. Because the machinery causing the shock wave 

moves, the shockwave rotates through a small angle A0 O . A second hologram is re- 

s 

corded. Because linear recording is used, the effect is the same as if the two holo- 
grams had been recorded simultaneously and the two refracted waves were propagating 
simultaneously from the vicinity of the shock wave . 

The figure exaggerates angles for clarity. The refracted waves occupy a spatial 
band that is easily determined from Snell's law of refraction. For the refraction from 
position 1, 


cos 0 

COS(0 g +0 X ) 


n r 


n, 


( 11 ) 


where the refractive indices on the consecutive sides of the shock wave are designated 
n^ and . For the refraction from position 2, 
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cos(0 s - A0 g ) 
cos(0 g - A0 g + e 2 ) 


( 12 ) 


n 

n 


2 

1 


Because the example refers to a gas, 0 ^ and 0 2 are small angles. By retaining 
only first-order terms in small angles in equations (11) and (12), the band occupied by 
the two refracted waves is given by 


2 Af 


ox 



^ s n 2 ~ n l 
A n 2 


(13) 


A large value for the refractive index change is given by 


n 2 - n^ = 0. 0003 

Also, n 2 is approximately 1.0. If the wavelength A of a helium-neon laser is chosen 
as 0. 633 micrometer, the spatial bandwidth of the equivalent transmittance of the dou- 
bly exposed shock wave is given by 


2 Af 


ox 


A0 . 


1/2 

mm • rad 


Clearly, a diffuser is needed to fill the entire numerical aperture of a television cam- 
era tube with a resolution of 18 per millimeter. 

Two additional simple calculations illustrate the other points made in this section. 
From reference 1, the depth of field 6 of the recording method is given by 


6 = 


1.03 

f 2 A 


(14) 


where f is the frequency resolution of the recording method. When A is 0. 633 mi- 
crometer and f is 18 per millimeter, the depth of field is 5 millimeters. A 2.5- 
centimeter-deep object could be resolved into about 6 longitudinal stations for true 
three-dimensional recording. 

Finally, suppose that we wish to use stereoscopic vision with an 18-per- 
millimeter recording. For an object point midway between the eyes, the minimum 
viewing distance S is given by 

S=— (15) 

2Af 
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where D is the pupil separation. When D is 6 centimeters, A is 0 . 633 micrometer, 
and f is 18 per millimeter, S is equal to 2 . 6 meters . Clearly, at that distance there 
would be little feeling of depth for a 2.5-centimeter-deep object. Methods for resolv- 
ing this problem are discussed in references 1 and 2 . 

Provided that the full resolution of a detector is available for recording the object 
wave, electronic recording of holograms is quite reasonable. A method for making 
the full resolution available is discussed in the next section. 


ELECTRONIC RECORDING OF HOLOGRAMS 

The main body of the analysis in this section is presented in references 1 and 2. 
The main purpose of this paper is to emulate in the spatial domain the time-domain 
processing discussed in these references. The processing technique can then be used 
with pulsed- laser holography. Hence, the time-domain analysis is repeated for quick 
reference. This section shows how to record a hologram that does not have the object- 
wave self-interference term and that also separates the virtual and real image-wave 
terms . Other sometimes-difficult-to-implement methods of electronic holography are 
described in references 12 to 14 and mentioned in references 1 and 2. 

An arrangement for electronically recording low-resolution holograms of semi- 
transparent objects is shown in figure 4 . A similar arrangement for recording holo- 
grams of reflecting objects is shown in figure 5 . This paper is concerned with semi- 
transparent objects of the type found in flow visualization; hence, figure 4 is primarily 
of interest. 

The key elements in figure 4 are 

(1) A beam splitter to combine the object and reference waves coaxially for low- 

resolution recording 

(2) An electro-optic phase modulator (ref. 15) for varying the phase of the refer- 

ence wave in time 

(3) A continuous-wave laser 

(4) An image dissector or any other scanned detector with an instantaneous re- 

sponse, as opposed to an integrating detector such as a vidicon 
The arrangement of figure 4 will implement the method known as heterodyne record- 
ing. In heterodyne recording, the reference wave contains different time-frequency 
components than does the object wave. Each component interferes with the object 
wave to generate a time-varying interference pattern. The pattern varies at a fre- 
quency equal to the difference between the frequencies of the object- and reference- 
wave components . An instantaneous detector is required to record the pattern at a 
position at a particular instant of time . 
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In heterodyne recording, the hologram must be scanned electronically . Scanning 
is discussed, for example, by Macovski (ref. 16). The object and reference waves 
will interfere to form a changing intensity pattern I(x, y,t). A scanning aperture 
S(x,y) is moved over this pattern at horizontal and vertical velocities V x and V . 
The instantaneous power from the scanning aperture is calculated by integrating 
I(x,y,t) weighted by the scanning aperture centered at its instantaneous position x c (t), 



dx dy S(x - x c ,y - y c )I(x,y,t) 


(16) 


The finite response of the system is represented by a detector impulse-response 
function h Q (t - t') that is nonzero only for t st'. The electrical signal upon scanning 
is given by 


f(t) = 



(17) 


Because of the need to scan back and forth in a raster, x Q (t) and y c (t) are not simple 
functions of time. However, on a scan line 


dx c (t) 

dt 


= V. 




> 


dt 



J 


(18) 


The finite duration of a scan line and the finite duration of a frame increase the 
signal bandwidth slightly in comparison with the contribution of I(x, y,t). In television 
practice, the frame is scanned in two fields to eliminate visual flicker. These effects 
add nothing to the analysis and are not considered further. 

The inverse process of writing a hologram pattern from a signal is similar. Sup- 
pose that the signal after processing, transmission, and reception is f^(t). A writing 
intensity is generated that is given by 


V) - 





t M )dt" 


(19) 
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The intensity is applied to a writing aperture W(x, y), where the writing aperture is 

scanned at horizontal and vertical velocities and Vy. At time t, the writing 

aperture is centered at x'(t),y'(t). The exposure at a point x, y is given by 

C/ ' c 


E(x,y) 



x c(t)»y - iw(t) dt 


( 20 ) 


The only other comment that will be made about writing a hologram in this manner is 

that, if V’ V v and V' ^ V , a scale transformation will result. 

^ ^ y y 

The detector response h Q (t) and the scanning aperture S(x, y) are chosen to be 
delta functions: 


h Q (t) = S(t) 


(21) 


S(x,y) = 6(x)6(y) 


( 22 ) 


This simplification is harmless if S(x,y) is narrow enough to resolve the desired ob- 
ject detail and if h Q (t) is narrow enough to pass a faithful replica of the scanned, time- 
varying intensity pattern. On a scan line beginning at 0,y Q 


x c (t) = V x t 


y e ( t )=y. 



(23) 


The signal is then given by 


f(t) = I(V x t,y 0 + V y t,t) 


(24) 


and is determined directly from the intensity distribution in the hologram pattern. 

Heterodyne recording is now easily explained. The object and reference waves 
are written in real notation as 


U(x, y) = 0(x, y) cos[cot - cp{x, y)] 
Ur( x , y) = s ( x >y) c °s[wt - r(x, y) + a(t) + 6 sin 



(25) 


where the circular frequency of the light wave is denoted by oj and a phase modulation 
term a(t) + 6 sin co^t has been introduced in the phase of the reference wave . Note 
that source retardation is negligible at the frequencies of interest. The reference wave 
can be phase modulated by using an electro-optic modulator. Temporarily, we equate 


I 
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a(t) to zero so that the remaining phase modulation 5 sin u>jjt depends on a depth-of- 
modulation parameter or modulation index 8 and a heterodyne circular frequency Wjj. 

By performing square-law detection to eliminate the dependence on co and by 
using the well-known Bessel function series (ref. 17) 




cos(6 sin Wjjt) = Jq(^) + 2 ^ J 2j £(< 5) cos 2 £w h 4 


1=1 




sin(6 sin Wjjt) = 2 2 J 2£+l( 6 ) sin [(2f + l)w H t] 


1=0 


(26) 


the hologram intensity is proportional to 


2 2 

I(x,y,t) = O (x, y) + s (x, y) + 2s(x,y)0(x,y) cos[<p(x,y) - r(x,y)] 


X 


J 0 (6) + 2 Yj J 2l< 6 > C ° S 2£C °H t 
1=1 


- 2s(x,y)0(x,y) sin[(p(x,y) - r(x,y)]i 2J 2i+i( 6 ) sin [( 2je + f (27) 


1=0 


J 


Note that the time variation contributed by scanning and the time variation con- 
tributed by phase modulation are quite independent. The depth-of-modulation parame- 
ter can be used to emphasize different harmonics. For example, by setting 5 equal 
to 2. 62, the first five orders of Bessel functions are 


J 0 (<5) = -0.106 
J L (6) = 0.465 
J 2 (6) = 0.461 
J 3 (S) = 0.239 
J 4 (6) = 0.0863 




> 


(28) 
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( 29 ) 


The first four terms in harmonics of o> R are 

I 0 (x,y,t) = 0 2 (x,y) + s 2 (x, y) - 0.212 s(x,y)0(x,y) eos[<p(x,y) - r(x,y)] 
^(x.y.t) = -1.86 s(x, y)0(x,y) sin(<p(x,y) - r(x,y)J sin co^t 
I 2 (x,y,t) = 1.84 s(x,y)0(x,y) cos[<p(x,y) - r(x,y)] cos 2co R t 
I 3 (x,y,t) = -0.956 s(x,y)0(x,y) sin[<p(x,y) - r(x,y)] sin 3u> H t 


Each of these terms contains a cross-interference term suitable for reconstructing 

the image waves . Because of the simple manner in which scanning affects the time 

dependence through equation (24), each spatial-frequency pair f , f contributes a 

x y 

corresponding time-frequency component 

h = V A + W < 30 > 


Generally, 


and 


^«1 


>1 

> 


^VxJ 


(31) 


In that case, each time-frequency component generated by scanning corresponds to a 
unique x-directed spatial-frequency component. 

By the proper choice of co R , each of the terms Iq, . . . can be placed in 
nonoverlapping bands and selected independently by band-pass filtering in the time 
domain. However, only Iq contains the self-interference terms; hence, the elimina- 
tion of Iq completes the heterodyne recording process and solves the overlap problem. 

If Af rx is negligible, the requirement for unqualified separation of all signals 
I 0 , Ip . . . is given by 


f„>3V Af 
H x ox 


( 32 ) 


For the equality sign, the bands for the first three signals (fig. 6) are 



I 


[*>• 2v x ^ox] 


[ 2V x A£ ox> ^x^ox] 

V 

rsv Af , 7V Af 1 
L x ox’ v x oxj 

If we assume that the scanning aperture does not resolve spatial frequencies 
greater than Af ox , the requirement on f^ can be relaxed somewhat . The separation 
requirement is 


‘H^x^ox 


(33) 


When the equality sign is selected, the bands are shown in figure 7 and are 


[o, V Af 1 
L * x ox J 

[ v x ^ox- 3V x ^ox] 

[ 3V x *ox" 5V x Afox] 


If fjj is decreased further, the bands overlap and the hologram constructed from such 
a signal yields overlapping multiple images. Equation (33) sets the lower limit on f^ 
in order to avoid overlap. 

The signals 1^ and Ig carry the complete hologram in either sideband. For the 
equality sign in equation (33), the lower sidebands are 


T- 



E ,V x Afox. 4V x Af 0 J 
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The corresponding signals are 


I x = 0.93 s(x,y)0(x,y) cos[jp(x,y) - r(x,y) - o> H f| 
I 2 = 0.92 s(x,y)0(x,y) cos[c?(x,y) - r(x,y) - 2o> H fj 


(34) 


The substitution x = V t has not been made explicitly. 

A hologram produced from either signal in equation (34) is an equivalent of an off- 
axis reference-beam hologram. For simplicity, choose the writing aperture W(x,y) 
and the writing response function h w (t) to be delta functions. If = V x and Vy = V , 
the hologram generated from 1^ contains the term 


I x (x,y) = 0.93 s(x, y)0 (x, y) cos 


<p(x,y) - r(x,y) - 


co H x 



(35) 


This hologram is equivalent to a hologram recorded with a reference beam that is off 
axis at an angle 6 given by 


sin 0 = — A s 2 Af A 

v ox 

x 


(36) 


The axis of the real- image wave appears at 0^, where 


2f 

sin0 ri = T7 iA - 4Af ox A 

v x 


(37) 


The virtual- and real-image waves are in the spatial-frequency bands [-Af QX , Af QX ] 

and [3 Af , 5 Af J , respectively. These results are summarized by figure 8. 
ux ux 

The following properties are needed to test the feasibility of heterodyne recording, 
fi 

If V„ Af^ v ft* 4x10 hertz, as with standard television (ref. 18), 

X ox 

f u ^ 2 V Af = 8x10® hertz 
H x ox 

For equality, 1^ is in the single sideband [4x10 Hz, 8x10 Hz] and is in the sin- 
gle sideband [12x10® Hz, 16x10® Hz] . Of course, the required frequencies and band- 
widths are reduced if V is reduced. The scanned system must have an h (t) narrow 
enough to respond to the maximum frequency . A television camera tube with the re- 
quired properties is the image dissector tube, which is essentially a photomultiplier 
in which an electron image is formed and scanned over a pinhole aperture. Finally, 
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to record the intensity given in equation (35), we still need a medium with resolution 
3 Af_,„. If the detector has a resolution of 18 per millimeter, the final recording me- 
dium should resolve 54 per millimeter. If resolution is a problem, fjj can be reduced 
to f^j after electronic recording so that 


f' 

sin 6 = — A > Af A 

v ux 

x 


2fl 


sine ri = 


H 


V 


A S 2 Af ox A 


The reconstructed waves are barely separated, but the required resolution is only 
36 per millimeter. 

Although phase modulation is an extremely convenient method for heterodyne re- 
cording of holograms and a good solution of the overlap problem, there is at least one 
other method of heterodyne recording in the literature (ref. 14). The unique feature 
of phase modulation is that it simplifies image processing. This property must also 
be emulated by a system that records pulsed- laser holograms. The use of continuous- 
wave electronic holography for image processing is discussed in the next section. 


IMAGE PROCESSING 

Methods of image processing as a means for extracting data from film are re- 
viewed in reference 19 . The theory of optical processing in continuous-wave elec- 
tronic holography requires a rather lengthy development, as presented in reference 1. 
A brief summary of the development appears in reference 2 . However, since we even- 
tually wish to show that the same operations can be effected in pulsed-laser holography, 
another brief summary is included here. 

Essentially, two operations are required in coherent optical processing: filtering 
of the spatial spectrum for image enhancement and cross correlation of the input wave 
for pattern recognition. Optically, both operations are performed by inserting a mask 
(spatial filter) in the spatial-frequency domain. In cross correlation the input wave is 
typically Fourier transformed by means of a lens (ref. 6) and then passed through a 
mask that has an amplitude transmittance proportional to the Fourier transform of the 
second wave. A second lens is used to invert the transform to complete the cross- 
correlation operation. This process is summarized in figure 9. 

The best-known application of the process shown in figure 9 is matched filtering. 
The mask is chosen to be the complex conjugate of the Fourier transform of the input 
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wave . Phase is canceled, and the output can be focused to a minimum-diameter spot 
by the second lens to show that the input is matched to the mask. For filter masks 
that can be expressed in the form 

H(x,y) = e-* a ( x ’ ^ (38) 

the entire process can be performed without lenses with heterodyne recording. The 
key operation is the multiplication of a hologram by an arbitrary phase distribution. 

If the phase modulation again contains a(t), the first three terms of the interfer- 
ence pattern are 

I Q = s 2 (x, y) + 0 2 (x, y) + 2J Q (6)s(x,y)0(x,y) cos[<p(x,y) - r(x,y) + a(t)f| 

I x - -4J 1 (6)s(x,y)0(x,y) sinf^x.y) - r(x,y) + a(t)] sin o> H t 
I 2 = 4J 2 (6)s(x,y)0(x,y) cos[<p(x,y) - r(x,y) + a(t)] cos 2o> H t 

The time t is associated with the position x„(t), y_(t) of the center of the scan- 
ning aperture. As in the description of heterodyne recording, a delta-function scanning 
aperture is chosen, and the scan line is to be horizontal. On a scan line, 

y = y G 

t=-i 

V x 

and a(t) is mapped into a function a(x,y). Hence, each cross-interference term in the 
hologram is multiplied by an aribtrary phase factor e ± J a ( x ’ y). 

However, the presence of a(t) in the phase increases the bandwidth with each sig- 
nal. Either the spatial-frequency content of the object wave must be decreased accord- 
ingly, or the scan velocity V must be reduced. The required restrictions are deter- 
mined as follows. The phase is to be controllable within the range - 7 r < a(t) < 7 r. The 
maximum frequency at which the phase is likely to be varied is V f , where f = 
Af_„ is the maximum x-directed spatial frequency in the object wave. Hence, the 

UX 

worst-case time variation of a(t) is given by 

a(t) = 7T sin STrV^^t (40) 

With the modulation index tt , at least the first four harmonics of V A f j atl must be 
counted on. Hence, a sideband must cover at least the range 5V x f xm . If the full res- 
olution of the detector is to be used, the scan velocity must be reduced by at least a 
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factor of 5. Interestingly, if the full resolution of the detector is used, the hologram 
finally constructed from the signal will represent an image wave of higher spatial fre- 
quency than could originally have been recorded. At this point, the system necessary 
to perform spatial filtering is available . 

In spatial filtering, an amplitude distribution 

, j<p(x 0 ,y 0 ) 

b ( x o’ y o) = °( x o’ y o> e ( 41 ) 


in an input plane is first Fourier transformed. A slight modification of the system 
shown in figure 4 will perform this transformation without lenses . The modified sys- 
tem is shown in figure 10. As derived in reference 1, the output signal is given by 



°( x 0 > y o> cos [,(V y o> ~~ (XX o + yy 0 > + a < t > " W H *] ^o dy o 


(42) 


where the correspondence between position and time is not shown explicitly. After it 
is constructed from the signal, the hologram is illuminated by an axially directed plane 
wave. The axes of the real- and virtual-image waves are directed at angles 
±sin -1 f H W x , respectively. 

The transformation must still be inverted by a lens. A better method is to modify 
the system depicted in figure 10 to use hologram- scale transformations. The scale of 
a hologram is easily changed in electronic recording by making the scan velocity 
V* , V* used to write the hologram different from the scan velocity V , V used to 

^ y ^ y 

read the hologram. If the hologram has a transmittance T(x,y) before the scale 
change, it will have a transmittance T(x/a,y/a) after the scale change, where the 
scale parameter a is given by 


a = 


V* = V£ 


(43) 


As shown in references 1 and 2, choosing an a of 0.707 makes it possible to use 
a diverging spherical wave to illuminate the object. The virtual image is placed at 
infinity and the real image is used as the output. The actual output for an input 
b( x 0 » y 0 ) is given by 

0.707 y\ (g) h*(0. 707 x, 0.707 y) (44) 


1^ = b* 1 0. 707 x - 


f H ks 
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where h is the inverse transform of the mask H . The output 1^ is the complex con- 
jugate of the desired output and appears magnified and off axis. However, none of 
these features is a problem in an application such as matched filtering. 

The recording and reconstruction processes are shown schematically in figure 11. 
The input transparency is placed at a distance S from the image dissector. The ref- 
erence wave originates at a distance 2S from the image dissector. The illumination 
of the subject transparency originates at a distance S in front of the transparency. 
Reconstruction is accomplished with a beam converging to a distance 2S in back of 
the hologram. The output appears at a distance S in back of the hologram but off 
axis. 

Lensless transformations and scaling can be performed with any electronic re- 
cording method. The key feature in using phase modulation with continuous -wave 
holography is that the phase 8 sin co-^t + a(t) can be inserted without reducing the 
available spatial resolution of the detector. The emulation of this method in pulsed- 
laser holography is the subject of the next section. 


ELECTRONIC RECORDING OF PULSED-LASER HOLOGRAMS 

In continuous-wave electronic holography, information is, in effect, stored in both 
the spatial domain and the time domain. As demonstrated, the entire spatial resolu- 
tion of the detector is then available for recording the object wave. In pulsed- laser 
holography, the time domain is no longer available. The only alternative is to try to 
control another attribute of the electric field. Optical holograms ordinarily record 
the magnitude and phase of a single polarization of the electric field. Hence, varying 
the polarization is another way to code information and is adopted as the approach 
here. The objective is to generate a signal like the signal represented by equation (39). 
Then, the overlap problem is solved, and an arbitrary phase modulation has been ap- 
plied. Furthermore, the entire spatial resolution of the detector must be available for 
resolving the object wave. 

To implement the system for pulsed-laser holography, the detector must have the 
following properties: 

(1) Because of the short duration of the pulse (20 nsec from a Q-switched ruby 
laser), the detector must be an integrating one. Most television camera tubes have 
this property. 

(2) The polarization sensitivity of the detector must be controllable, one resolution 
element at a time. That is, given two orthogonal directions on the resolution element, 
the sensitivities to intensity of light polarized in those two directions must be indepen- 
dently controllable. 
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Unlike the continuous -wave case, where the image dissector has the desired proper- 
ties, there seems to be no detector that has a polarization sensitivity controllable on 
an element-by-element basis. There are physical phenomena, however, that might be 
exploited to realize this property. These phenomena include variable dichroism in 
crystals (ref. 20) and magnetic domain control (ref. 21). 

An experimental arrangement would consist of modifying a few items in figure 10. 
The required modifications are as follows: 

(1) The image dissector is replaced with a detector of the type previously defined. 

(2) The reference wave is circularly polarized by inserting a quarter-wave plate 
in the reference-wave path. (The laser emits a linearly polarized beam. ) 

(3) The object wave is linearly polarized at equal angles with respect to the con- 
trol directions. 

(4) The phase modulator is replaced by a ''black box" for controlling polarization 
sensitivity. 

Maintaining a linear polarization of the object wave is relatively easy in flow- 
visualization applications. A diffuser can be constructed that does not depolarize the 
illuminating beam, and the weakly refracting object does not disturb the polarization. 

The mathematical proof of this system is presented here . The circularly polar- 
ized reference wave has the form 

U r = s x cos(wt - r) i + Sy sin(wt - r) j (45) 

A A 

where i and j are unit vectors along the x- and y-directions, respectively. The 
vector U r and the scalars s x , Sy, and r are possibly functions of position. Note 
that s x = Syj however, subscripts are retained to identify the polarization. A real 
notation will be used throughout this section; hence, oj is the circular frequency of 
the light. The linearly polarized object wave is given by 

U Q = o x cos(u>t - cp) i + Oy cos(u>t - cp) j (46) 

Again, o x = Oy although subscripts are retained to identify polarization. The vector 
U Q and the scalars o x , Oy, and <p are possibly functions of position. The square- 
law detectable intensity I is given by 

I = I f («o + u? + 2 u 0 u> 

Integral over many 
optical cycles 


2 2 

o + s 
x x 


+ s x° x cos(<p-r) + 


2 2 
o + s 

— ^ ^ + s o sin(<p - r 

9 y y 


(47) 


where t is the integration time . 
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A major requirement is that the polarization sensitivity not increase the resolution 
requirement of the detector: The entire spatial resolution is to be available for resolv- 
ing the object wave . To check that this requirement is met, let the detector consist of 
an array of square apertures . Weight each aperture with a polarization sensitivity 
sin [a(X, Y)/2] in the x-direction and a polarization sensitivity cos [a(X, Y)/2] in the 
y-direction, where the coordinate of the center of the aperture is denoted by X, Y. 

The aperture is assumed to have a side length L. The signal read from a resolution 
element during a scan is proportional to the expression 



The subscripts have been dropped from s and o. Here, the aperture must be small 
enough to reproduce the maximum -frequency component in sin [a(X, Y)/2] by sam- 
pling. 

The aperture must also be small enough to resolve the maximum-frequency com- 
ponent in so cos(<p - r). 




dy so cos(cp - r) « so cos (cp - r) x y 


(49) 

L 


where the subscript X, Y denotes evaluation of the function at the center of the aper- 
ture . The point to note is that the entire resolution of the detector is available for re- 
solving the object wave. The resolution conditions apply separately to the two factors 
in the cross-interference terms of equation (48). A detector to resolve the pattern 


. 2 
sin 



so cos {cp - r) 


directly would need a resolution equal to the sum of the maximum frequencies in the 
two factors. 
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Noting that the resolution requirement is satisfied and keeping that requirement in 
mind, for simplicity assume a delta -function aperture . The signal read during a scan 
is then proportional to 

S = o 2 + s 2 + sin 2 so cos (cp - r) + cos 2 j ^ x? y -^j so sin(<p - r) (50) 


The use of squared trigonometric functions has ensured that the self-interference 
2 2 

terms o + s are not multiplied by a sensitivity factor. The signal can be written as 

S = o 2 + s 2 + — cos (cp - r) + — sin(<p - r) 

2 2 


SO so 

- cos a(x,y) — cos (cp - r) + cos a(x,y) — sin ((p - r) (51) 
2 2 


The last two terms in equation (51) will yield reconstructed waves containing an arbi- 
trary phase in the manner of equation (39). We show that the phase can be used to 
place these terms in bands that do not overlap the self-interference terms. 

Select a(x,y) = <5 sin 27rfx, where f is the maximum spatial frequency that can be 
sampled by the aperture. As discussed, the frequency f will be the maximum spatial 
frequency that can be resolved in the object wave. For an x-directed scan velocity 
V , the corresponding time signal from the cross-interference terms containing an 

A 

arbitrary phase factor is 

S = -cos(<5 sin 27rfV x t)s(V x t,y)o(V x t,y)[cos <P(V x t,y) - r(V x t,y)Tj 

+ cos(<5 sin 27rfV x t)s(V x t,y)o(V x t,y) sin [^(V x t,y) - r(V x t,yj] (52) 


After the usual expansion, the first three terms are 

S 0 = -J 0 (S)s(V x t,y)o(V x t,y) [cos(<p - r) - sin {cp - r)J 

5 1 = -2J 2 (S)s(V x t,y)o(V x t,y) c°s(4 5 rfV x t) [cos(^ - r) - sin(<p 

5 2 = -2J 4 (6)s(V x t,y)o(V x t,y) cos(87rfV x t) [cos(^ - r) - sin(^ 




- r)] 

- r)] > 


(53) 
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In keeping with the term "heterodyne, " we define 

“H - 2rtV x 

We also recognize the correspondence between time and position by writing V t as x. 

A 

The terms become 


S 0 = -J 0 (<5)s(x,y)o(x,y) [cos(<p - r) - sin(<p - r)] 

5 1 = -2J 2 (<5)s(x,y)o(x,y) [cos(<p - r) - sin(^ - r)] cos 2co H t L 

5 2 = -2J 4 (<5)s(x,y)o(x,y) [cos(<p - r) - sin (<p - r)J cos 4co H t 


(54) 


These terms are similar to those in equation (27) or (29) with the following exceptions: 

(1) Odd multiples of are missing. 

(2) The cos (cp - r) and sin(<p - r) are both present, although this feature does not 
affect the spectrum or the reconstruction. 

The various terms of the spectrum occupy the following bands (fig. 12): 

(1) Self-interference terms occupy the band (0, 2fv ). However, the detector 
does not resolve spatial frequencies greater than f so that the actual signal corre- 
sponding to these terms is in the band (0, fV ) . 

(2) The Sq term occupies the band (0, fV ) . 

(3) The S n terms occupy the bands 

(fVx, 3fV x ) 


<3fV x , 5fV x ) 



Clearly, the terms multiplied by cos 2No>jjt, where N is an integer, do not overlap 
the self-interference terms or cross-interference terms in the time domain and can 
be selected by filtering. 

At this point, the pulsed-laser system has been proven. It is interesting to note 
that the criterion for bands not overlapping is that 
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(55) 


H 


Af V 
ox x 


However, if f = Af QX is the maximum resolvable spatial frequency, the equality sign 
must apply: 


f H = M ox V x 


(56) 


In continuous-wave holography, the criterion according to equation (33) was 


f H- 


2 Af V 
ox x 


(57) 


That criterion cannot be satisfied in the pulsed- laser case. The missing odd multiples 
of cojj make electronic recording of pulsed- laser holograms feasible. To implement 
image processing, the arbitrary phase distribution a(x, y) is applied by varying the 
polarization sensitivity . 


CONCLUDING REMARKS 

For the benefit of the cursory reader who might be interested in recording holo- 
grams electronically, a few concluding remarks are in order. These conclusions are 
easily drawn from the theory or the examples. 

Electronic recording of holograms is intended primarily as an information- 
processing method, not as a viewing method. For a single recording, the width of a 
nearby object is limited to the width of the detector. Even then, the eyes must be at a 
large distance from the hologram for stereoscopic viewing. However, the development 
of a display for three-dimensional viewing is not entirely precluded. There are elec- 
tronic plotters that will display orders-of-magnitude more resolution elements than 
can be resolved by a detector at one time. This resource can be used to make viewing 
easier in the following ways: A wide-angle view can be recorded as a series of views. 
An array of detectors can be used to record an extended object. To implement stereo- 
scopic viewing, an object wave can be multiplied by a random-phase diffuser during 
recording, or the image-wave terms can be separated from the self-interference 
terms by filtering, stored in a computer memory, and then multiplied by the random- 
phase diffuser. 

The development of a system for electronically recording a hologram really con- 
sists of several projects: The feasibility of a detection method must be demonstrated. 
An appropriate detector must be developed or found in the pulsed- laser case . The 
value of the detection method in image processing must be proven. Attempts to write 
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a hologram from a signal result in another set of projects. However, the variety of 
materials for storing holograms and the higher resolution attainable in writing a holo- 
gram from a signal should make these projects easier. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, November 30, 1978, 

505-04. 
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APPENDIX - SYMBOLS 


a 

a(t) 
a(x, y) 

b ( x 0 »y 0 > 


D 

d 

E 


f H 


ox 

Af 


ox 


rx 

^rx 

f t 

f(t) 

V) 

H(x, y) 
h 


V*> 

v*> 

I(x,y,z),I(x,y,t) 


3 

a 

j 


scale parameter 
arbitrary-phase modulation 

arbitrary phase to be added to phase of object wave 

amplitude distribution at an input plane of an image-processing 
system 

pupil separation 

distance from recording plane to local origin of surface 

exposure 

spatial frequency 

heterodyne frequency 

central spatial frequency of object wave 

half -bandwidth of object wave 

central spatial frequency of reference wave 

half-bandwidth of reference wave 

time frequency 

electrical signal 

signal following processing, transmission, and reception 
amplitude transmittance of filter 

inverse transform of filter mask (impulse response function of an 
optical system) 

detector impulse response function 

impulse response function for writing a hologram 

intensity 

writing intensity 

unit vector in x-direction 
Bessel function of first kind, n n order 
imaginary designator 
unit vector in y- direction 
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k 


L 

n l’ n 2 

°( f x’ f y> 

°x 9 °y 

o(x,y) 

P i 

R( f x> f y) 

r(x,y) 

S 

S(x,y) 

s x’ s y 
s(x, y) 

T(x,y,z) 

3T 


t 

U(x,y, z) 
U R (x,y,z) 


V x’ V y 

W(x,y) 

x,y,z 

x c ( t ),y c (t) 

y 

6 

<5(x), 6(y), <5(t) 
e 


wave number 

length of square scanning aperture 
refractive index 

Fourier transform of object-wave amplitude 

magnitudes of x- and y-components of object wave 

magnitude of object wave 

instantaneous power from scanning aperture 

Fourier transform of reference-wave amplitude 

phase of reference wave 

signal or distance from hologram 

scanning aperture 

magnitudes of x- and y-components of reference wave 

magnitude of reference wave 

amplitude transmittance 

Fourier transform 

time 

phasor representing an object-wave amplitude 
phasor representing an image-wave amplitude 
real vector amplitude of object wave 
real vector amplitude of reference wave 
scan velocities in x- and y-directions 
writing aperture 
Cartesian coordinates 

instantaneous position of center of scanning aperture 
proportionality constant 

either depth-of -modulation parameter or depth of field 
Dirac delta function 

angle of reconstructed virtual-image wave relative to optical axis 
angle of reconstructed real-image wave relative to optical axis 
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6 „ angle of a shock wave 

b 

A 9 angle of rotation of a shock wave 

0p02 angles of refraction 

A. wavelength 

<P ox angle of object-wave axis relative to x-direction 

(pjtx angle of real- image axis 

<P rx angle of reference-wave axis relative to x-direction 

<P X angle of propagation of plane wave relative to x-axis 

(Py angle of propagation of plane wave relative to y-axis 

<p(x,y) phase of object wave 

co circular frequency of light 

cOjj heterodyne circular frequency 

Superscripts: 

* complex conjugate 

' equivalent 
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(a) Schematic of hologram-recording process in x,z-plane. 
(Reference wave has negligible spread. ) 


Real-image 

wave 



Figure 1. - Schematics of hologram-recording and reconstruction processes. 


-Cross-interference 




-fv 


(a) Recorded with off-axis reference wave. 

Figure 2. - One-dimensional spatial-frequency spectra of hologram. 


(b) Recorded with off-axis reference wave when single side- 
band will be retained. 


32 


Position 1 



Figure 3. - Double exposure of refraction from rotating shock wave. 



Figure 4. - Arrangement for recording low-resolution hologram of trans- 
parency. 
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Figure 5. - Arrangement for recording low-resolution hologram of 
solid object. 
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Figure 6. - Spectrum of nonoverlapping bands for unqualified separation of signals for f H 
“V- Af 
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Figure 7. - Spectrum of nonoverlapping bands when frequencies exceeding maxi 
mum object-wave frequency are not resolved and f H * 2V X Af ox . 


Beams as viewed head-on at 
large distance from hologram 



Figure 8. - Reconstruction of hologram recorded by heterodyne method. 
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Figure 9. - Spatial -filtering operation using pair of lenses. 
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Figure 10. - Spatia (-filtering operation using ienseless transform and heterodyne 
processing. 
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